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We investigate the excitation function of directed flow, which can provide a clear sig- 
nature of the creation of the QGP and demonstrate that the minimum of the directed 
flow does not correspond to the softest point of the EoS for isentropic expansion. A novel 
technique measuring the compactness is introduced to determine the QGP transition in 
relativistic-heavy ion collisions: The QGP transition will lead to higher compression and 
therefore to higher compactness of the source in coordinate space. This effect can be 
observed by pion interferometry. We propose to measure the compactness of the source 
in the appropriate principal axis frame of the compactness tensor in coordinate space. 

1. Motivation 

The primary goal for the investigation of heavy-ion collisions is to test the equation 
of state (EoS) of hot and dense matter far off the ground state, especially with view on 
possible phase transitions, e.g. to the Quark-Gluon-Plasma (QGP) [||]. Indeed, collective 
flow phenomena are sensitive indicators for thermodynamically abnormal matter 0. In 
the case of a first-order phase transition to a QGP, an isentropic expansion proceeds 
through a stage of phase coexistence which should lead to signatures in the observables. 
The first ideas to investigate this phenomenon occured in the mid-seventies [0. In this 
paper we investigate the excitation function of directed fiow, as well as the compactness 
in heavy-ion collisions. 



2. Model 

To investigate quantitatively the experimental observables, we perform 1-fiuid and 3- 
fluid (3-1-1 )-dimensional relativistic hydrodynamic calculations. That is, we solve numer- 
ically the continuity equations for the energy-momentum tensor, df^T'^'^ = 0, and the net 
baryon current, d^Nj^ = 0. Detailed discussions of (3+l)-d numerical solutions for hy- 
drodynamical compression and expansion can be found e.g. in [^. We shall employ two 
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different equations of state for P{e,p): 

i) A relativistic mean field (RMF) hadron fluid corresponding to baryons and an- 
tibaryons interacting via exchange of massive scalar and vector bosons, plus free 
thermal pions; the parameters of the Lagrangian are fitted to the ground state of 
infinite nuclear matter, in particular the nuclear saturation density, the energy per 
particle, and the incompressibility. 

ii) The same EoS as in i) for the low density phase, but supplemented by a Bag Model 
EoS with a bag constant B^^^ = 235 MeV for the quark-gluon (QGP) phase. The 
phase coexistence region corresponding to this first-order transition is constructed 
employing the Gibbs' condition of phase equilibrium, Phg(T, ^ib) = Pqg(T, fis), 
where T and fiB denote the temperature and the baryon-chemical potential, re- 
spectively. For example, for p = we find Tc ~ 170 MeV, while at T = phase 
coexistence sets in at p ~ 4.6po- A more detailed discussion of these EoS can be 
found in 0. 

Further, we employ the three-fluid model with a dynamical local unification procedure [0]. 
The three-fluid model treats the nucleons of the projectile and target nuclei as two different 
fluids, since they populate different rapidity regions in the beginning of the reaction. The 
same holds for the newly produced particles around midrapidity, which are therefore 
collected in the third fluid. Thus, the three-fluid model accounts for the non-equilibrium 
situation during the compression stage of heavy-ion collisions. The coupling between the 
projectile and target fluids is calculated assuming free binary A^A^-collisions 

The unification of fluids i and j consists of adding their energy-momentum tensors and 
net-baryon currents in the respective cells. 



and common values for e, P, p and u'^ are obtained from T^^jgg^ = {e + P) u^u'^ — P g^" , 
^unified — pu'^, and the given EoS P = P{e,p). The local criterion for unification is 
{Pi + Pj)/P > 90%. Here, Pij denotes the pressure in T-^j , and P the pressure in T^j^jgg^. 

3. Directed flow and softest point of the EoS 

In order to measure the EoS, i.e. the pressure P(e,p) as a function of energy density e 
and baryon density p in the local rest frame of a fluid element, the transverse momentum 
in the reaction plane, p^, is investigated. This quantity is proportional to the pressure 
created in the hot and dense collision zone [0]: 



The pressure P is exerted on a transverse area element A_i . Directed flow has therefore 
been proposed as a measure for the pressure and a possible "softening" of the EoS @J^. 

Fig. |l| shows the excitation function of directed flow p'^^ /N calculated in the three- 
fluid model in comparison to that obtained in a one- fluid calculation 0. The one fluid 
calculations show that for increasing bombarding energy, the flow, ~ px, first increases 



Tr{x)+T^\x)=TZU^) 



(1) 




(2) 
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(for a collision without phase transition), as the compression and thus the pressure grow. 
At large E^^^ the time span of the collision decreases, diminishing the flow again. The 
flow is thus maximized at some intermediate bombarding energy. In the case with a phase 
transition the decrease of the flow is much more rapid than for the purely hadronic fluid. 
The reason for this is not that Cg, the isentropic velocity of sound, vanishes but rather 
geometry: The compactness and tilt-angle B are different in the calculation with phase 
transition, and this leads to a different initial condition for the subsequent expansion (see 
below and ||10|| ). After passing through a local minimum at E^^^ ~ 5A GeV, the directed 
in-plane momentum reaches a second local maximum around E^^^ ~ 10 - 20A GeV. This 
is the point where the compressed matter first becomes hot enough (over a large volume) 
to "respond" with small pressure gradients. 

Due to non-equilibrium effects in the early stage of the reaction, which delay the build- 
up of transverse pressure |TI|, the flow in the three- fluid model is reduced as compared 



to the one-fluid calculation in the AGS energy range. Furthermore, the minimum in 
the excitation function of the directed flow shifts to higher energies. The case without 
dynamical unification yields the least amount of stopping and energy deposition, while 
the one-fluid calculation has instantaneous full stopping and maximum energy deposition. 
The three-fluid model with dynamical unification lies between these two limits; it accounts 
for the limited stopping power of nuclear matter in the early stages of the collision and 
mutual equilibration of the different fluids in the later stages. Most importantly, the 
three-fluid calculations predict an increase of p'^"/N towards ^'Lab — GeV, if indeed 
a phase coexistence with small Cs occurs. Data at that energy has recently been taken, 
and should prove very useful to pin down the onset (or absence) of a first-order phase 
transition in the AGS-SPS energy domain. 

First order phase transitions "soften" the EoS i.e. P increases slower with e and p 

than in the case without phase transition. This corresponds to a reduction of the isentropic 
speed of sound, Cg, as compared to that in the interacting hadron fluid. However, as shown 
in Fig. ^, this happens only if the entropy per net baryon in the central region is large. 



i.e. if the ratio T/^b is not too small |]T0|. In the three-fluid model that is due mainly to 
the kinetic equilibration of the decelerating baryon dense projectile/target fluids with the 
midrapidity fluid of secondary particles, which leads to considerably larger s/p than in 
one-fluid hydrodynamics. However, at the energy corresponding to the minimum of the 
directed flow the specific entropy is rather small, s/p < 10; i.e. in the AGS energy domain 
matter is rather baryon dense but not very hot. Consequently, the EoS is not soft (i.e. 
the isentropic velocity of sound is not small), even if mixed phase matter does occur. 

Fig. ^ shows the time-like component of the net baryon current in momentum space 
{Px — Piong plane) for Pb-|-Pb-collisions at 6 = 3 fm. One clearly observes the directed 
in-plane flow (before the collision, there is no matter at Px ^ 0). However, in the left panel 
there is almost no momentum of baryons in the upper left or bottom right quadrants, 
where p^ ■ Piong < (except for two "jets", see below). The central region passed through 
the phase coexistence region at high s/p, with a rather small average Cg, and isentropic 
expansion of the highly excited matter is inhibited. Note the difference to the expansion 



pattern observed for lower energies, right panel of Fig. |^, where (c^) is not small fTO 



The slope of (px/N) at midrapidity is shown in Fig. ^ as a function of beam energy. 
Experimental Data are shown as well One observes a steady decrease of Ey = 
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d{px/N)/dy up to about top BNL-AGS energy, where the flow around midrapidity even 
becomes negative due to preferred expansion towards ■ piong < 0. The "overshoot" 
towards negative slope is due to the small incompressibility in the top AGS energy region, 
and the rather early fluid uniflcation employed here. However, such a behavior can not be 
observed in the data. A less steep decrease of Fy could be achieved in the three-fluid model 
by a more stringent uniflcation criterion (i.e. later uniflcation) or early kinetic decoupling 
on the hadronization hypersurface. At higher energy, Eiab — 40^4 GeV, we encounter 
the expansion pattern depicted in the left panel of Fig. |: flow towards p^ ■ Piong < 
can not build up! Consequently, Fy increases rapidly towards Eiab = 20 — 40^4 GeV, 
decreasing again at even higher energy because of the more forward-backward peaked 
kinematics. Note that the increase of the slope is due to the absence of the "anti-flow", 
see Fig. ^. In any case. Fig. ^ shows that it will be difficult to see the effect of the possible 
phase transition in Fy. The double-differential in-plane cross section. Fig. ^ appears more 
useful. 



4. Compactness 

Measurement of the compactness is a promising new tool to observe the onset of the 
phase transition. It relies on measuring the shape of the source, which is uniquely related 
to the pressure and density of the system in the compression and expansion stage of 
the nucleus-nucleus collision. The compactness can be identified via interferometry: The 
illuminiation of the baryon source by the pion radiation is subject to experimental scrutiny 
via pion interferometry Fig. ^ illustrates the basic idea. It shows the baryon density in 
the reaction plane for the EoS without (i) and with (ii) phase transition, respectively. One 
clearly observes the higher compression in the case with phase transition. As indicated 
above, the onset of the transition to quark matter at a given incident energy Ef^J^ leads 
to higher compression p/po than for the case without transition. Now, as pV ~ nR^Lp 
must equal 2A by virtue of baryon number conservation, the longitudinal thickness L of 
the compressed matter is approximately proportional to 1/p. Thus, a transition to quark 
matter leads to a more compact system, just as quark matter stars are more compact than 



pure neutron stars |jT5|. Of course, in heavy- ion collisions that expectation is based on 
the behavior of relativistic compression shocks rather than hydrostatic and gravitational 
equilibrium. 

In particular, we study the compactness in the reactions Au-|-Au at impact parame- 
ter 6 = 3 fm for the energy -E^^T — 8^ GeV. The compactness is defined as the ratio 
of the smallest to the biggest in-plane eigenvalue of the configuration space sphericity 
tensor, which we define as the second moment of the net baryon current. On fixed-time 
hypersurfaces we have 



F- ■ 



J e (p(f ) - Pent) . (3) 



We apply an additional density cut p > pcut in the integral to discard spectator matter. 
In the future the cuts and the hypersurface will have to be adapted to the experimental 
conditions. However, this is not crucial for understanding the effect. 

The three eigenvalues fn are the solutions of the cubic equation det(Fjj — fSij) = 0, and 
the eigenvectors e„ follow from solving the linear systems of equations {Fij — fnSij)el^ = 0. 
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In terms of the eigenvalues /„ and orthogonal eigenvectors of F, the compactness tensor 
can be written as F = fiCi ® el + /2e2 62 + /sCs ® 63. In diagonal form, F specifies 
an ellipsoid in configuration space with principal axis along e„ and radii y/f^. Cigar-like 
patterns, oriented along the z-axis, would lead to /i > f2 = fs, el = e^, (^2 = (^x, 63 = e^. 
On the other hand, a "pancake" / "lensil" shape corresponds to /i < /2 = /s- The tilt 
angle is determined from the scalar product of (the longitudinal direction in the lab 
frame) with the vector e„ corresponding to the biggest of the eigenvalues /„. As already 
indicated in the introduction, we find very different eigenvalues for the two equations of 
state. The calculation with transition to quark matter corresponds to higher compactness 
of the baryon distribution. That is, the compactness tensor is much flatter (nearly a factor 
of two !) in the model with (ii) than in the model without (i) phase transition. Moreover, 
our (3+l)-dimensional expansion solutions show that after the compression stage the ratio 

of the in-plane radii \J~f2lT1 remains much smaller in the case with phase transition, cf. 
Fig. ^. Note also that the extremely fiat shape of the baryon distribution means very 
small curvature of the surface, which in turn will lead to a strongly "bundled" emission 
of hadrons from the (almost planar) rarefaction wave or deflagration shock by which the 



dense droplet decays; see also the discussion in 

The eigenvalues of the compactness tensor allow to measure directly the density increase 
in the high density stage of the reaction, if a phase transition occurs. Care must be 
taken that the impact parameter range investigated constitutes a moderately small bin 
of centrality values. One should keep in mind that the compression factor is affected 
by the incompressibility dP/de evaluated on the shock adiabat (in the one- fluid model), 
not along a path of fixed specific entropy. Therefore, the incompressibility needs not be 
equal to the isentropic speed of sound. The latter is not much reduced in the presence 



of the phase transition at AGS energies ITQ], because of the high net baryon density and 
relatively low temperature. 

Lisa et al. have recently proposed a new interferometry analysis, which could be used 
to observe this change in the compactness directly. The developed method is quite robust 
and incorporates other interesting information as the configuration space tilt angle, which 
nicely complement the momentum-space flow angles. It avoids cuts in tilted ellipsoids, 
which are not analysed in the appropriate rotated frame, and where the excentricity and 
the RMS-radii are much less distinct for the two different equations of state. 
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Figure 1. The excitation function of directed flow Figure 2. Pressure P as a 

p'^^ /N for Au + Au collisions at impact parameter function of energy density 

6 = 3 fm. Dotted lines (open symbols) are results e for different entropy per 

from one-fluid dynamics; triangles are for a purely net baryon ratios for the 

hadronic EoS, circles are for an EoS with phase tran- EoS with phase transition 

sition. Solid lines are calculated with the three-fluid to QGP. 
model, with (large circles) or without (small circles) 
dynamical unification. All three-fluid calculations are 
performed with phase transition. 




Figure 3. Net-baryon density in momentum space. Pb(40 AGeV)+Pb (left; s/p ~ 20) 
and Pb(8 AGeV)+Pb (right; s/p < 10) at 6 = 3 fm. 
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Figure 4. The slope of the directed in-plane momentum per nucleon at midrapidity for 
Au-Au-coUisions at b=3 fm (three- fluid model with dynamical unification), the experi- 
mental data shown is from ||13|| . 

8 AGeV, b=3fm, 1 -Fluid Limit no PT 8 AGeV, b=3fmi, 1 -Fluid Limit with PT 
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Figure 5. Baryon density in the reaction plane for i^L^b — GeV, at time tcm = 6 fm/c. 

Left: EoS without phase transition. Right: EoS with phase transition. 



8 AGeV, b=3fm. 1 -Fluid Limit 
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Figure 6. Ratio of the in-plane radii y/2//i for the RMF-EoS without phase transition 
and for the case with transition to quark matter. 
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